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ABSTRACT 

We investigate the coupling efficiency of parametric downconversion light (PDC) into single and multi-mode 
optical fibers as a function of the pump beam diameter, crystal length and walk-off. We outline two different 
theoretical models for the preparation and collection of either single-mode or multi-mode PDC light (defined 
by, for instance, multi-mode fibers or apertures, corresponding to bucket detection). Moreover, we define the 
mode-matching collection efficiency, important for realizing a single-photon source based on PDC output into a 
well-defined single spatial mode. We also define a multimode collection efficiency that is useful for single-photon 
detector calibration applications. 

1. INTRODUCTION 

Earliest studies of parametric downconversion (PDC) addressed problems in fundamental physics, while more 
recent studies target applications such as quantum metrology^ and quantum information.^''^ While both of 
these areas make use of two-photon light, they are distinct applications that present different requirements for 
that light. The stringent requirements of these applications are driving researchers to optimize the PDC process. 
For these efforts to succeed, a clear theoretical framework is needed. 

PDC produces a quantum state of light with a two-photon field description. However, if only one photon 
of the pair is measured, the source exhibits purely thermal statistical behavior, given by intrinsic multi-mode 
photon production. We can however introduce a certain degree of coherence. By measuring one of the photons, 
we prepare the other photon in a specific state. The prepared state will be pure only if we project the first 
PDC photon (called also the heralding photon) into a pure state. In each of the above PDC applications, we 
prepare one photon by measuring its twin. Thus for optimization of the process, it is crucial to have a proper 
definition and measurement of the efficiency of that preparation and of the related mode-matching. We present 
two different models to define and optimize the two-photon-mode preparation and mode-matching efficiencies. 
They are distinguished by how the heralding photon is collected. One uses a multimode spatial filter or a bucket 
detection system, while the second uses a single-mode fiber. A bucket detector is a multimode detector where all 
the modes are detected jointly. Hence the information about the location of the detected photon (or equivalently 
in what mode the photons were detected) is " erased" . We obtain different dependencies of the efficiency on the 
pump parameters in these two arrangements. This is particularly important for two specific applications: the 
calibration of a single-photon detector and the realization of a single-photon on demand source (SPOD).*'^ 

2. BIPHOTON FIELD 

We consider a two-photon wavefunction, written as^ 

IV') = / d^pid^p2diidi2^'(pi,P2,ti,i2)|lp2,t2>|lpi,ti)> (1) 

where px^i represents the transverse positions of the two photons at the instant 1^.2 and $(^1,^2,^1:^2) is 
the biphoton field. Its calculation is analytic only for first-order approximation of the transverse wavevectors 
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(assuming the pump, signal, and idler have narrow transverse angular distributions we can adopt the paraxial 
approximation), and is obtained by performing the Fourier transform of the pump angular distribution and 
the phase-matching function with respect to the pump transverse k-components and the signal k-components. 
Perfect transverse phase-matching is also assumed. The rcsiilt derived in icfJ is 
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where t = ti — t2 and IIdl{t) = 1 for < t < DL and elsewhere. The subscripts s, i, p indicate the 
signal, idler, and pump. Oi^^ are the central emission angles (in a small angle non-coUinear approximation), and 



i,pi^i,s,p,(t^)^i,s,p/c describe the directions of the central intensities of the wavevectors. The terms. 
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p — — ciliui — — account for the effects on the refractive indexes ( np_s,i{i-^p.s.i, 4') 

expanded around the central frequencies {fls,i), and around the phase-matching angle <po) of the pump and the 
signal due to the pump angular spread, which are responsible for a small deviation from the phase-matching 
angle (po- D = - ^"'(^')"'/^ |a, + d»sK,0)c^s/c |^^ ^j^^ differential phase velocity between the signal and idler 
photons in the crystal. D = and Afg = for type I degenerate phase-matching. The pump-beam transverse 
field distribution is assumed to be Gaussian with a waist of Wp at the crystal. We also assume that the pump 
propagates with negligible diffraction inside the crystal. In the following, the idler and signal space-time positions 
are indicated by 1 and 2, respectively. 



3. SINGLE MODE EFFICIENCY 

Envisioning the optics setup as the unfolded scheme of Fig. 1,^ the source is described by the propagation of a 
coherent mode defined by ipimiPz) through an optical element with impulse response function hi{p3, pi), through 
the non-linear crystal where the mode gets transformed according to the phase matching function $(pi, p2, ^i, ^2), 
and collected eventually by hs{p2, Pi)- The actual collected mode will then be given by the field ^pim{pi)- 

The coincidences measured at the positions 3 and 4 are then C34, the overlap between the PDC field and 
both of the preparation or collection modes, while the single coimts C3 and C4 measure the overlap individually 
between the biphoton field and each of the preparing or collecting modes: 



■-34 



-I 



d^Pld'^P2d'^p3d'^pA'^{pl,tl,p2,t2)hi{pi,P3)hs{p2,p4.)^l'rn'{p4.)'Plrn{p3) 



I 



dtidt2d Pi 



d^Pld^P2d^P3 ^{pi,ti, P2,t2)hi{pi, P3)hs{p2, PiMmiPs) 



Ca = 



d<idt2d^P3 d^pid^p2d^/!?3 ^{pi,ti, P2,t2)hi(pi, P3)hs{p2, p4)^*lm{Pi) 



(3) 



Here we distinguish between the single-mode matching efficiency calculated by 

C34 



Xm = 



and the single-mode preparation efficiency 
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Figure 1. Unfolded picture of an object corresponding to the idler channel of a PDC source, generated in a nonlinear 
crystal NLC of length L. The object (on the left) is a coherent source when the idler beam is prepared by single-mode 
fiber. 

Eq. (4) is the appropriate efSciency to optimize when the final goal is a single photon source on demand emitted in 
defined photon modes, while Eq.(5) is the efficiency measured when single photon detector calibration is required. 
Specifically in the paper, we calculate the efficiency in a perfect imaging configuration, namely hs{p2. Pi) = 
5{p2 — M4P4) and hi(p3, pi) = d{pi — M^p^) (lenses have infinite aperture with magnification Mj+2). The lenses 
are arranged to place the preparation and collection beam waists, = 250 //m at the crystal, with guided 

Gaussian field modes ifiQ{pj+2) — v fi^^ ^^P — ^ with j ~ 1, 2. The spatial coherence of the single 

guided modes in the signal and idler arms should ultimately match the overall spatial coherence of the two-photon 
states. By explicit calculation the single-mode matching efficiency is then given by 
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and B ^ {Np ~ Ns + KpOsfwl-^ + {-Np + Ns + KpOifwl^ + K^{di + Osfwl- In the thin crystal limit this 
reduces to 



as first calculated in Refs..^'^° The single mode preparation efficiency is given by 
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Also in this case we can calculate easily the thin crystal limit 
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Note that 



XM = ^|xfxf. (12) 

In this case the information of the location of the prepared and detected modes is included in the measurement. 
In Fig. 2(a), we plot the theoretical prediction of Xp for ^ 5 mm crystal of LilOa versus the collecting waist 
Wo, 2 for fixed preparing waist Wo,i = 250 fim at the crystal, for different pump waist configurations. As it 
appears there is an optimum collecting waist and the best efficiency is obtained for the larger pump waist, xm 
behaves as xp, except that the highest efficiency is obtained when the preparing and collecting waist are the 
same (wo, 1=^0,2)- This can be clearly observed in Fig.2(b)where the optimum collecting waist is plotted versus 
the pump waist for different prepared waist Wo,i- The optimum Wo,2 (giving the highest efficiency) approaches 
Wo,\ as the pump waist increases, and this approach is faster in the xm case, as expected. 

4. MULTI-MODE OR BUCKET DETECTION EFFICIENCY 

In the case of multi-mode or bucket detection efficiency the schematic is Fig. 3, where we overlap the PDC wave 
function with the impulse response function of the optical system, and eventually we average over the spatial 
distribution of multi-mode fibers or spatial filters (apertures) or bucket detector, erasing the information about 
the spatial location where the photons where collected. More specifically, in Fig. 3, the source is considered 
incoherent with finite transverse distribution 73(^3). The collecting modes are then represented by the spatial 
filter and described by ^4(^4). The coincidences measured at the positions 3 and 4 are then C34, and C3 and C4 
are the single counts: 
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The 2-photon multi-mode matching efficiency is then defined by^^ 

C34 

and the 2-photon multi-mode preparation efficiency is defined by 
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Figure 2. (a) Plot of xp versus the collecting mode waist for various pump waist Wp — 150, 200, 400, 600 nm at Wo.i = 250 
/xm (b) Plot of the optimum waist Wo,2 versus the pump waist for various preparing waists Wo,i in the case of xm (dashed 
line) and xp (solid line). 



We calculate the multi-mode matching and preparation efficiencies, assuming the incoherent source is com- 
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pletely determined by the functions 7^4.2(^^+2) = e , where the impulse response function is the Dirac 

delta as before. Calculating the multi-mode matching efficiency we obtain: 
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In the thin crystal limit rjM becomes 
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Also in this case the multi mode preparation efficiency is given by 
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Figure 3. Unfolded picture of an object corresponding to the idler channel of a PDC source, generated in a nonlinear 
crystal NLC of length L. The object (on the left) is an incoherent source prepared by a multi-mode fiber, and the image 
is collected by a multimode fiber. 




Figure 4. Plots of rjM (dashed) and xm (solid) for Wo,i — Wo.2 = Wo (solid lines) and various pump waists, Wp — 
150, 200, 400, 600 ^tm, versus Wo- 



with thin crystal limit 
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In the multi-mode approach presented here, the preparation and collection modes can be thought of as 
spatially filtering or selecting the multi-mode input light. As Fig. 4 demonstrates, the predictions made by this 
model yield different results than the single-mode model. The multi-mode model predicts that, for a fixed pump 
waist, the maximum multi mode-matching efficiency is obtained when the fiber-defined collection mode (at the 




Figure 5. Plot of rjp^M versus the collecting mode waist for various pump waist, Wp = 150, 200, 400, 600/im and for fixed 
preparing waist Wo.i = 250 ^m. 

crystal) is large, i.e. all the pumped crystal volume is in a region of unit collection efficiency of the spatial filter 
system. With the single-mode model, one of the fibers acts as a source of a single-mode beam that propagates 
through a spatial filter (in this case the pumped crystal volume) to the other fiber. The maximum mode-matching 
efficiency is achieved with a large pump waist, with respect to the preparation and collection beam waist at the 
crystal. If the pump waist is smaller than the fiber-defined collection beam waist, the mode-matching efficiency is 
reduced. In this case, the efficiency is maximized with Wo.i = Wo.2 — but this presents the practical difficulty 
of having, and aligning, exactly the same selected modes. Here, the differences between the two models are more 
evident. In fact, when the collection/preparation waist is much greater than the pump waist, r]M asymptotically 
goes to 1, while xm goes to 0; in the opposite condition {wp ^ Wo), rjM — + ^ and xm 1- 

In Fig. 5, we plot the theoretical prediction oirjp^M for a 5 mm crystal of LilOa versus the collecting waist Wo,2 
for fixed preparing waist Wo,i = 250 /im at the crystal with a range of pump waists. As it appears for r\p there 
is not an optimum collecting waist and the highest efficiency is reached faster for the smaller pump waist. The 
T/M presents instead a finite optimum waist for smaller pump waist, but the highest efficiency is never reached, 
because modes are matched only on average, given the statistical nature of the multimode model. 
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Figure 6. Unfolded scheme for the preparation of the heralding channel in a single mode, while photons are collected in 
a bucket detection scheme by an aperture or a multi-mode fiber placed in the focus of a lens 

5. SINGLE MODE PREPARATION AND MULTIMODE COLLECTION 

EFFICIENCY 

When we prepare the source according the unfolded scheme of Fig. 6 in a single mode, while we collect in a 
multi-mode, the single counts are C3 in Eq.(3), and the coincidences are: 
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The single-mode preparation and multi-mode collection efhciency is 



eP = (21) 

We explicitly calculate it with a mode preparation given by a single-mode fiber in a perfect imaging configuration 
(lens with infinite aperture and magnification M3 arranged to place the collection beam waist, Wo 1 at the crystal), 

described as a Gaussian field, (pio{p3), while collecting in a multi-mode with a Gaussian aperture 74(^4) = e 
The signal optical system is hs{p2,P4) = exp[—ik{p2 — P4)^)], the free space propagation, with k = n/{Xsd) and 
d the propagation distance. The single mode preparation and multi-mode collection efficiency is given by 
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Figure 7. Plot of ep versus the collecting iris aperture for various pump waist, Wp — 150, 200, 400, 600/im 

B' = {JVp-Afs + KpOsfwl^wlk^ + {-Np+Ns + KpO.fiwl^ +wl + Pw^wl^wp + K^kV,^^w*{9, + O^f, and 
C" = {w1 1 + Wp)^ + k'^w'l iWplw'^Wp + wl iiD^ + vol ivo^). In the thin crystal limit this reduces to 

We observed that in this case ep depends slightly on L. We plot in Fig. 7, the theoretical prediction of ep for 
a 5 mm crystal of LilOs versus the iris aperture w for fixed preparing waist Wo,i = 250/im at the crystal, for 
different pump waist configurations. In this case we have an intermediate behavior between xp ryp, the 
highest efficiency is obtained faster for larger pump waist. We underline that the comparison with the following 
experimental results show a qualitative analogy with this curve, because the k parameter plays an essential rule 
in the estimate of the efficiency. Specifically the quantitative value of k in the simulation does not correspond 
to the experimental values. 

6. PRELIMINARY EXPERIMENTAL RESULTS 

We give here some preliminary experimental results and compare them with theoretical predictions, showing the 
agreement. The simulations and measurements are done exciting PDC in a type I phase-matching configuration 
with a pump wavelength of 351 nm and by a LilOa crystal length of 5 mm. Pump waists at the crystal were 
ranging from from 150 /im to 600 /im. 

Some measurements were done with the trigger arm (heralding channel) coupled with a single-mode fiber. 
The lens images the minimum waist at the crystal. The signal arm was either coupled with a single-mode fiber in 




Figure 8. Plots of ep data versus the iris diameter (a) for a range of pump waist with fixed Wo = 250 and (b) for a 
range of pump waists with Wp ^ Wo- (a) shows measurements performed at lEN, where the efficiency is uncorrected for 
the detector dead-time, while (b) shows measurements performed at NIST, corrected for dead-time. The solid line are 
fits obtained with the expected theoretical curve in the thin-crystal limit reported in Eq.(24). 



a perfect imaging configuration or with a multimode-fiber placed in the focus of the coupling lens, where an iris 
at the lens selects the collecting modes. Other measurements were performed in a similar setup (6 mm LilOa) at 
the National Institute of Standards and Technology (NIST) with the trigger beam coupled with a single-mode 
fiber and the signal coupled to the detector by a single lens, with an iris selecting the modes. In particular 
two experimental configurations were chosen. In one configuration the trigger mode was fixed to maximize the 
singles rate, while in the other configuration the heralding channel waist was matched to the pump waist at the 
crystal. 

In Fig. 8 we show the experimental values of the single-mode preparation and multi-mode collection efficiency 
(ep) versus the collecting iris diameter, for the two experimental configurations, i.e. for fixed Wo = 250 /im and 
(b) for fixed Wp ~ Wo and a range of pump waists as indicated. In both cases the maximum efficiency is obtained 
for larger pump waists as predicted by the theory in Fig. 7. 

We point out that the theoretical prediction of the multi-mode collection and multi mode preparation ef- 
ficiency (j]p) versus the equivalent Gaussian filter collection aperture Wo,2 would instead give the maximum 
efficiency for smaller pump waists. However, while the single-mode preparation and matching efficiency (xm,p) 
versus the collecting waist Wo.2 present a maximum for properly matching the three modes, 77p,j\/ may not reach 
the highest efficiency even in this case. 

Fig. 9 shows the experimental single counts from a heralding channel coupled with a single-mode fiber versus 
the fiber waist at the crystal, for various pump waists. Data are fitted by theoretical curves that for single 

mode coupling are C-i — , — , — - — ; here the data clearly matches the single- mode propagation 

because we have an optimum preparation waist that maximizes singles rate (note that the maximum singles rate 
are for a smaller pump waist). In a multi- mode configuration the single rates would just increase with the iris or 
filter aperture. In Fig. 10 we plot the experimental data of the uncorrected efficiency in the case of single mode 
fiber placed on both detection arms for a 5 mm LilOa versus pump waist at the crystal. Measurements were 
performed for two different preparing and collecting mode waists configuration, that is the apparent waists at 
the crystal were Wo,i = HO fJ.m and Wo,2 = 50 fim and reversed. The data are fitted by the Eq.(ll) and Eq.(12), 
in the thin crystal configuration, this is because of the very small size of the waists, where the theory for long 
crystal is no longer reliable. The agreement shows the correctness of model proposed. 




Figure 9. Plots of single counts C3 versus Wo for a range of pump waists. The maximum counts are for the smallest 
pump waist. 




Figure 10. Plots of the single mode preparation efficiency xp s^nd matching efficiency xm versus the pump waist Wp for 
two preparing and collecting mode configurations, Wo,i = 111 fim and Wo.2 ~ 50 /im the upper curve and Wo,i = 50 fim 
and Wo, 2 = 111 t-ini the lower curve. The curve in the middle is xm for "^0,1 = 155 /im and 'Wo,2 = 70 fi. The solid lines 
are fit obtained by the proposed theoretical model. 



7. CONCLUSIONS 



In summary, wc have presented an analytic model to quantify the mode preparation and matching efficiency 
in terms of adjustable experimental parameters with the goal of optimizing single-mode collection from PDC 
sources. In addition, we have presented an alternative scheme that may have more validity for multi-mode 
collection arrangements. We reported preliminary experimental results, supporting the validity of the proposed 
multi-mode collection model. 
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